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In this third part of our calculation of the QCD NLO corrections to the 
photon impact factor we combine our previous results for the real corrections 
with the singular pieces of the virtual corrections and present finite analytic 
expressions for the quark-antiquark-gluon intermediate state inside the photon 
impact factor. We begin with a list of the infrared singular pieces of the 
virtual correction, obtained in the first step of our program. We then list the 
complete results for the real corrections (longitudinal and transverse photon 
polarization). In the next step we define, for the real corrections, the collinear 
and soft singular regions and calculate their contributions to the impact factor. 
We then subtract the contribution due to the central region. Finally, we 
combine the real corrections with the singular pieces of the virtual corrections 
and obtain our finite results. 

I. INTRODUCTION 

This paper represents the third part of our investigation of the photon impact factor in 
next-to-leading order QCD. First results of the calculation of NLO corrections to the photon 
impact factor have been published in [l|] and 0. In [||] we have computed the QCD NLO 
corrections to the process 7*+^' — >■ {<l<i) + (l which lead to the virtual corrections of the photon 
impact factor. The second paper, 0, contains the process 7* + g ^ {(I^IQ) + 5', leading to 
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the real corrections of the photon impact factor. Results are given for the helicity-summed 
squared matrix elements, with the longitudinally polarized photon in the initial state. In 
the present paper we will complete this part of our calculation by listing also the results for 
the transverse polarization of the photon. 

The main purpose of the present paper is the combination of singular pieces of virtual 
and real corrections. To begin with the former ones, in we have listed the results for 
the virtual corrections which contain both the finite and the singular pieces. As the first 
step we therefore have to separate the sum of all singular 1/e terms from the finite part. 
In the next step we turn to the real corrections. Infrared singularities will come from those 
kinematic regions of the qqg system where the gluon is either soft or coUinear with the quark 
or with the antiquark. We compute the soft and the collinear approximations of our matrix 
elements, and, by subtracting these divergent contributions from the full matrix elements, 
we form finite combinations. Finally, by combining the singular pieces of the real corrections 
with those coming from the virtual corrections, we obtain further finite terms which have 
the form of the Born approximation and can be computed analytically. As the main result 
of this paper, we present finite expressions for the qqg contributions to the photon impact 
factor. As the fourth and remaining part of our program we will be left with the task of 
evaluating numerically the finite integrals obtained in this paper, as well as the finite pieces 
from the virtual corrections. 

This paper will be organized as follows. After the definition of the impact factor (sec- 
tion H) we first (section HD list the infrared singular pieces of the virtual corrections. We 
then turn to the real corrections and present the complete results for the real corrections, 
both for the longitudinal and for the transverse photon (section [TV] ). In the following sec- 
tion we compute the collinear and the soft limits of the real corrections. We discuss in 
detail the definition of the impact factor correction in section |VJ, where the subtraction of 
the central region, i.e. the leading log s term, is discussed. In section |VII| we complete our 
program, by combining the singular pieces of the real corrections with those of the virtual 
corrections and by demonstrating their infrared finiteness. In a final section we summarize 
our results and give a brief outline of the remaining part of our program. 



II. DEFINITION OF THE IMPACT FACTOR 

We consider the elastic scattering of two particles A and B with momenta pa and in the 
Regge limit s — oo, where the momentum transfer t = q^ is kept fixed. In leading order in 
as we write the scattering amplitude Tab in the form: 

J-AB\s,t - q } -IS j ^27r)o-2^ r2 (g - r)2^ ■ 

The impact factors and for the external particles (quarks, gluons, or photons) are 
of order a^. The exchanged gluons carry transverse momenta r and (q — r) respectively. 
Taking the cut amplitude we define the total cross section for the scattering of A and B via 
the optical theorem as 

= ilmT^»i(., < = 0) = / (g^tS'^g' . (2) 
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Here we have dropped the tilde symbol on the impact factors to denote the additional s- 
channel cut. This equation can be used as the definition of the leading order impact factors 
^^^^ and Note that we have summed over the colour indices of the t-channel gluons 

imphcitly. The total cross section can also be written as 

^aI = YsI \^AB\'d<^^ (3) 

where l/2s denotes the fiux factor (in the high energy limit), A^^^ the matrix element for 
the scattering process A + B A' + B' with the exchange of a gluon, and d(j) stands for the 
phase space of the final state A' + B'. For our case of interest, 7*g (or 7*7*) scattering. A' 
denotes a qq pair, B' a quark (or a qq pair). 



A 



B 




FIG. 1. Kinematics for A + B^A' + g + B' with initial state momenta pa,Pb and final state 
momenta k\,i, fe^. An S^^-reggeon is exchanged, emitting a real gluon with momentum £ which is 
assumed to be separated from the other final state particles by large rapidity gaps. 



In the high energy limit, the scattering amplitude Mas for the scattering A + B A' + B' 
(with the invariant masses of the particles or systems A' and B' being finite) (Fig. |^ without 
the extra produced gluon) is described by the exchange of a reggeized gluon and can be 
written in the following form: 

g 

Mab = ^^A-,A' 

Here, uj(t) stands for the Regge trajectory of the gluon, and F^^^, and ^b^b' 
particle-particle-reggeon vertices. The index a denotes the colour of the exchanged reggeized 
gluon. Eq.(^) exhibits the factorization property of Regge theory: when written in complex 
angular plane, the residue of the gluon Regge pole can be written as a product of two vertex 



+ 



— S 



^ B^B' 



(4) 
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functions, one for the incoming particle A A', the other one for B B' . In the total 
cross section formula this factorization property leads to the factorization of impact 
factors. Since (|^) is valid not only in leading order, also the validity of the impact factor 
representation is rather general. However, when going beyond the leading order, we have 
to replace the t-channel propagators of the exchanged gluons by reggeon propagators. In 
leading order, we put uj{t) in @ equal to zero and obtain 

A/f(0) _ -nW.a 2S (o),a , . 

■'^^AB — ^ A->A'~^ B->B' ■ K"^) 

Inserting this expression into (^, we arrive at (|^). 

In the next order, a^, eq.(]^) receives several corrections. In order to illustrate the general 
pattern it is instructive to start from (|^). There are three new contributions on the right- 
hand-side of @. The first one comes from the same intermediate state A', B' as in the 
leading order case, but there are higher order corrections inside the vertex functions Ta->a' 
or Tb^b' of the scattering amplitude A + B A' + B' (or its complex conjugate). A 
second correction is due to the reggeization of the exchanged gluon: from it follows 
that in next-to-leading order reggeization provides an extra term in Mab proportional to 
r^l^'^/a;(t) ln(s/ — t)r^^J^^,. These two corrections will be referred to as 'virtual' corrections 
to the impact factor. The third correction comes from the production of an extra gluon in 
the intermediate state ('real' corrections): A + B^A'-\-g + B' (Fig. |l|). In the high energy 
limit, we divide this intermediate state into configurations with one or two large rapidity 
gaps: in the former case, the extra gluon belongs to the fragmentation region of A or B, in 
the latter case to the central (or: multiperipheral) region. 

Correspondingly, we separate the integral over the phase space d(f). Begin with the 
configuration with one single rapidity gap. If Na' [Nb') particles in the final state belong 
to the system A' {B') and carry momenta k\, (/k'^/), we may write the phase space in D 
dimensions explicitly as 

d(P = rf0^,rf0B'(27r)^5(^)(pA +PB- E^k\' - E.k'B') (6) 

with 



d k At 



^A' 



and 



JB' 



Introducing the explicit integration over the reggeon momentum r with a Sudakov decom- 
position r = a^q' + l3rP + via one extra delta function and 

d^r s dar dPr d^^'^r . . 

[9) 



(27r)^ 2 27r 27r (27r) 



D-2 
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we may write 



sdPr J, sdar 1 



■>A'— ^VB 



27r 2tt 2s(27r)^- 



(10) 



where d(j)A' = #A'(2vr)^5(^)(pA + r - E.^a') and = (27r)^rf0B,(5(^)(j9B - r - ^ .A;^,) 
are the usual phase space measures for the final state particles belonging to the systems A' 
and B' respectively. For an intermediate state with two rapidity gaps (gluon in the central 
region) we have two t-channel reggeon momenta, r and r'. Denoting the momentum of the 
produced gluon hj i = r — r' = agq' + Pep + i± eq. ([T0|) will be generalized to: 



(11) 



1 2 sdPr sdar sdPr' sda^' d^ '^r d^ ^r' 

where 

= (^^^^(^')(2vr)^5(^nr -r' + i) (12) 

is the phase space measure for the gluon in the central region. 

Returning to the impact factor representation (|^), the NLO corrections due to the extra 
gluon in the s-channel can be grouped into two parts. The gluon in the fragmentation region 
of A or leads to corrections of the impact factor $^ or <I>b, respectively, whereas the gluon 
in the central region represents the first rung of the BFKL ladder (in the leading-logarithmic 
approximation) 0. Together with the other NLO corrections mentioned before, the entire 
sum of NLO corrections to the total cross section can be cast into the following form: 



+/(c)^^'^^^(^/^')2"'^(^')>^' 

+ l-(^/^o) / (^^^$r^/C.ea,(r,r')^$2^ (13) 

Here /Creai(^, ''"') denotes the square of the gluon production vertex, the BFKL kernel (some- 
times referred to as the 'real' part of the BFKL kernel, i.e. it does not contain the gluon 
reggeization). Later on it will be convenient to rewrite the reggeization of the gluon (third 
term on the right-hand-side): putting ln(s/r^) = ln(s/so)+ln(so/T^) we combine the In (s/sq) 
piece with the real part of the BFKL kernel: in the sum of both terms the infrared singularity 
in the limit r — r' — > drops out. Note the factorization of eq. (|l^): as indicated after (H), 
this feature is a general consequence of Regge theory and (with suitable generalizations of 
the i-channel gluon propagators) is expected to hold in arbitrary order perturbation theory. 
For the impact factor we state the general definition: 



=yE / {n-^A'^'A^A^dc^A^-^ (14) 
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The average symbol () stands for the sum over color and helicity in the intermediate state 
A'. In the following we will use the short hand notation 

^A = J2f \^A^A'\'d<j>A''-^, (15) 
A' ^ 

In particular, in the square of the matrix elements we have included the colour projector 

In this paper we are interested in the NLO corrections to the photon impact factor. 
The strategy follows from the discussion outlined in this section. The NLO corrections 
can be divided into the qq and qqg intermediate states (named virtual and real corrections, 
respectively). In the former case we need the NLO corrections to the 7* qq vertex, r^V_^gg. 
They contain infrared divergences (in dimensional regularization: poles in e) and are listed in 
An additional (infrared divergent) contribution comes from the gluon trajectory function 
in the third term on the right-hand-side of (plsf) . For the qqg contribution we need the vertex 
r^*^qqg, in the Born approximation: since, by definition, the gluon lies in the fragmentation 
region of the virtual photon, we have to divide the phase space of the produced gluon. In 
section |V1| we present a detailed discussion of this rather subtle issue: we will show that the 
separation of the fragmentation region introduces the scale Sq, in the last term of ([T3|). The 
fragmentation region contains soft and collinear divergences. When combining them with 
the infrared divergences of the qq state (including those from the gluon trajectory function), 
the infrared divergences cancel @] . It is the main purpose of this paper to find a separation of 
the fragmentation from the central region which allows to define infrared finite combinations 
of virtual and real corrections to the impact factor. 



III. SINGULAR TERMS OF THE VIRTUAL CORRECTIONS 

For future reference, and in order to define our notation we briefly summarize the Born level 
impact factor and list the divergent terms of the virtual corrections. The amplitude T^*q for 
the scattering •y* + q' ^ qq + q' (Fig. ^ has been calculated in to next-to-leading order 
in ttg. We write T^*q as an expansion in the strong coupling g 

Tr, = g'TiZ + g'T\l . (16) 
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The momenta are labelled as in Fig. ^. The kinematical variables we use are the centre- 
of-mass energy s of the incoming virtual photon and the incoming quark, the virtuality of 
the photon, Q^, the Bjorken variable x = Q'^/2p-q, the momentum transfers t = r^, = k"^, 
tf, = {q — k — r)"^, and the invariant mass of the outgoing gg-system. In addition we 
use a Sudakov decomposition of the momentum k with respect to the light cone momenta 
q' = q — xp and p with 2p- q' = s , k = aq' + l3p + k±. Two dimensional transverse momenta 
are denoted as = —k^. 

The Born impact factor is obtained from T^?g. To leading order we may write 

^(0) ^ ■p(0),a 2s (0)^^ 

■J-^q i-^'^qq ^ ^ qq \^ ' ) 



with the Born level vertices 



Tf),a = ].u{p-r,\q,) /AMp, A,), (18) 



where 



H^ = u{k + r,\) ^ }i ^X'^v^q-k.X') , (20) 
= u{k + r,\) ^{4- /(f- /) ^ X'v^q - k, A') (21) 

and A" are the generators of the colour group SU{Nc). Squaring the Born level vertices 
(p!8|) and averaging (summing) over incoming (outgoing) colours and helicities we obtain for 
the coupling to the reggeized gluons with colours (a, b) , contracted with the colour singlet 
projector (cf. eqs. ( p^ and (|T5|)): 



|rJ?P = ^^^- (22) 



The 7* gg-vertex |ry.^_^g^p = \T^^}^g^{k, a)p for a longitudinally polarized 7* coupling to 
a quark flavour / reads 

|r?L Ji ^ .VW^le^aHl - afQ^ (5^ - ^) ' . (23) 

where D{k) = k"^ + a{l — a)Q'^. Here the exchanged reggeized gluons must be in the colour 
singlet. For the sum over the two transverse photon polarizations we have 



|r!,».L,,|?. = 2VNr^le^a(l - a) [a^ + (1 - a)^ (^ - 



(24) 



The LO impact factor follows from the definition in eq. (plSj), from the two-particle phase 
space measure 

sdl3r da d^~'^k 



2n 2a{l - a) (27r)^-i 



(25) 
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and from the Born level vertices ( ]23| , Expanding in the strong coupling, 



we write the Born level impact factor as 

d^-^k da 



$(0) 



,(0) 



The last equation defines the integrand of the LO impact factor: 



J2;L(a, k- r, Q) = ^^^-^^ a\l - afQ^ 



(27r 



1 



D{k) D{k + r) 



(26) 



d''-'kdal2-T,L{(^^k;r,Q). (27) 



(28) 



l2-^T{o!,k;r,Q) 



(27r 



D{k) D{k + r] 



(29) 



Because of the symmetry of the squared matrix element under the exchange q ^ q, namely 
a ^ 1 — a, k ^ —r — k, we have the relation X2(a, k; r, Q) = X2(l — a, —r — k; r, Q). 



Turning to the NLO part of (|l^), we split T^*q into a finite part F and a divergent part 
D as 



7^(1) 

7*9 



F -u ^(0) ^ 

7*9 "I" //1„\2-e 7*9 ' 



(471 



cr = 



r( i + 6)r^(i-6) 
r(i-26) 



(30) 



The divergences are given in dimensional regularization and appear as poles when the space- 
time dimension D = A — 2t approaches D = A. Here we explicitly include the divergent 
contribution from self-energy insertions on all four external quark lines 



D 



SE 



-2Cr 



(31) 



labelling poles from ultraviolet divergences explicitly by euv- As usual, the SU{Nc) invari- 
ants are denoted by Ca = Nc and Cp = {N^ - l)/{2Nc). 

The divergent part can then be written as the sum of three terms 



D — Dj*^gq + + Dgq , 



(32) 



following the expansion of the Regge ansatz (^) in powers of a^. The three terms denote the 
divergent parts of the vertex r^i-*,^^, of the Regge trajectory of the gluon, u'^^^t), which is 



accompanied by the large leading logarithms in s, and of the vertex T'^gq , respectively. They 
have the form: 



2Cf 



3Cf Po , 2Cf t ( j^2\ , (t a{l — a)t ^ , sq 



+ ^ + ^ln(-M2) + 
euv e 



In- 



M2 



+ In- 



(33) 



D,. 



Ca 
t 



■5 , -s 
m h m 

So So 



(34) 
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and 



2Cp 3Cp /Sq 



qq 



. ^ Ca , So 

-t) + — In — . 

e —t 



(35) 



Here, we have separated infrared and ultraviolet divergences, making use of the formulae 
given in Note that, compared to the Regge ansatz (H), we have changed the energy scale. 
In (0), we have introduced the energy scale sq: instead of lns/(— t) we write lns/(— t) = 
In s/ So + In so/{—t) and absorb the second term into D^^gg and Dgg in (|33| , ^). Let us first 
consider the infrared singular terms. For the virtual NLO corrections to the photon impact 
factor we start from 



(1, virtual) 

7* 



s. 



ab 



da 



.{l),a p{0),b* 



a) 



+99 T 



*qq""rqq 



sdPr 
2^ 



+ c.c. . 



(36) 



Making use of eqs. ( ^Ol) and (^3|) — the latter one without the UV pole which we will discuss 
later — we can immediately deduce the IR divergent part of the virtual correction to the 
impact factor: 



(1, virtual) 



IRdivergent 



cr 



(47r) 



dk / dal2{a,k;r,Q) 



Cr 



[An] 



2-e 



— r21ogso + 21oga(l-a)-21ogM21 + — 
e -' e 



-- + 41ogM^ 
e 
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(37) 



Below we will show that these divergent terms cancel against the real corrections. 

Turning to the ultraviolet divergent pieces in D^^gg and Dgg, they are proportional to 
/3o = (llA^c ~ 2?2j)/6, and in |^ it has been shown that they are removed by the renormal- 
ization of the strong coupling. It is instructive to repeat the argument in more detail. The 
ultraviolet part of (|33D , including constants, is obtained by adding, in eqs. (56), (57), 
(61), (62), (67), plus one half of eq. (60): 



4 

9 Cr 



.(0),a 



(An) 



UN, 1 

Uf — 

6 3 ^ 



1 



2iV. 



(38) 



Except for the last term which is due to the wave function renormalization of the massless 
outgoing quarks, all ultraviolet poles are proportional to (r^)"*^. It will be convenient to 
slightly modify the last term: instead of using, for the renormalization of the outgoing 
massless quarks, the decomposition (|3l|) , we write 



SE 



-2Cf 



(r^) 



2\-e 



(39) 



This allows to simplify 



g^cr (o),a (r 

(47r)2- -y'-"^ euv 



-(3o 



(40) 
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When addressing, in the next step of our program, the calculation of the finite terms of 
the virtual corrections, we will have to remember that the change (|3^) in the wave function 
renormalization introduces a new finite term 

^\v^}'Z,,C,Hr'). (41) 



(4vr) 

In order to remove the UV divergences we have to renormalize, i.e. we add the Born 
approximation with the replacement: 



(42) 



with the renormalization constants taken in the MS scheme. When both contributions are 
added, the ultraviolet poles cancel, and we are left with 

/r»;^l=-^n»>_«Aln;^. (43) 



Note that all dependence upon the renormalization scale /i is contained in the logarithm 
which introduces the scale dependence of the strong coupling constant, agir^y Correspond- 
ingly, in the photon impact factor the /i-dependence is contained in the contribution: 



2j<.,uv,^_«,^,„, (44) 



as it is required by the renormalization group equation. Our derivation might suggest that 
is the 'natural' momentum scale of in the photon impact factor: however, it is clear 
that we can write 

.2^(l;UV) _ _ ^(0) ( n.'^ L n 1. 



^^^r" = -i*^\/5oln^+/3oln^J (45) 

i.e. by simply redefining the /^-independent part of the impact factor we can switch from 
one momentum scale of a^, a.s{r^^^ to another scale, ^^(Q^). 

Stating that ( ^41) is the only //-dependent term in the impact factor implies that the 
/f-dependence drops out of all the remaining contributions to the impact factor. In fact, 
we silently dropped a /t already in eq. (^I]), anticipating the scale invariance of the finite 
terms which are left after the cancellation of the infrared divergencies of real and virtual 
corrections. This can be understood as follows. Without having it written explicitly in [|T|, 
the infrared divergences in dimensional regularization appear as poles in e, accompanied 
by the power of some dimensionful scale, Sj. In a one- loop calculation the sum of infrared 
divergent terms from real (R) and virtual (V) corrections generally reads 



Adding real and virtual corrections, the coefficients Ai and Bi are sums as well, 

Ai = AY + , Bi = BY + Bf. (47) 
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Expanding (|^) in e we obtain 

TiR 



— In Si + — In Si H 5^ In s^. 

e 2 e 



(48) 



Introducing the scale n with the renormalized coupling as a factor n'^'" and reexpanding 
results in 

In In^ /i^ ' 



IR 



+ 



A. 



1 - 2 1/11 1 2 

- + ln/i I In Si + —In Si 



+ - + ln/i2 Si - Si In Si. 



(49) 

Further below (section |V11|) we will show that in the sum of virtual and real corrections all 
e-poles cancel: this means that Ai = and Bi = 0. Hence, we see that the cancellation 
of infrared divergences implies that also the In /i^ terms cancel. This allows us to write e.g. 



the logarithm in (|4l|) with a dimensionful argument. We did the same with the results of 
diagram 14 in after expanding divergent terms of the form s^^. 



IV. REAL CORRECTIONS 



r 



q — k — £ q 




Q p p — r ^ 

FIG. 3. Kinematics of the process 7* + g' — > qqg + q'. 



We now turn to the real corrections. We study the process 7* + 0' — ^ {(1(19) + Q with transverse 
photon polarization (Fig.|^); the case of the longitudinal photon has been studied in 0, and 
we follow the procedure outlined in this earlier paper. The diagrams are shown in Fig. ^. We 
use the Sudakov decomposition of the gluon momentum i = a^q' + jd^p + £j_ and introduce 
the following abbreviations : 

ai = a , (50) 
0.1 = {1 — a — ai) , (51) 
a2 = {l-a), (52) 
0:2 = (a + a^). (53) 
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FIG. 4. Feynman diagrams for 7* + <?' ^ qqg + q'- 



The propagator denominators, that occur in the calculation, have the following form: 
Di = {k + i + r-qf 



\2 "1 -,2 



at 



= - as - — (fc + r)^ 

= {k + ^ + rf 



0L\ OLi 



1 



£'4 = (A; - qf 



1 



{at{k + £) + ai£f 



at 



^ -aiQ^-^(k + ry-^{k + e)'-(e-r) 



a\ 



ai 



Da = k' 



c,,Q^-e-^{k + £f-^£^ 



ai 



D^^{k + ef 



ai 



ai 



{k + £f. 



(54) 



(55) 



(56) 



(57) 
(58) 
(59) 

(60) 
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We present the results for |A4gqc,p, the squared matrix element, averaged (summed) over 
the helicities and colors of the incoming (outgoing) quarks and summed over equal photon 
polarizations; thereby we made use of X]i=j=i 2 (^t ^-L)(fr ^J.) = ^-^^ Our results can be 
simplified by expressing them in terms of the matrix elements for longitudinal photons; for 
this reason we also recall the results for longitudinal photon polarization 0. The matrix 
elements for the process 7*g —>■ {qqg)q can be written in the following factorized form: 



I A// - |2 = I 



^"i xab |p(0) 1 2 



(61) 



and by using eq.(p2D we extract the squared vertex |r!^°''_^gggp (summed over helicities and 
color in the qqg state). Following the notation used in we list our results in correspondence 
with the Feynman diagrams: 



= 4 e\)g' E (^4 + ^^^^ + ^4 + ^4) ' (62) 



(63) 



c 



aiai(l -e) {k + £ + rY {k + r)' 



(64) 



+ (^{k + rf + ai Q') {2 (k + r){k + £ + r) + ai {{k + rf + Q')) ^ 



+ 



2{k + rf . ,i 
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ai {k + r)^ (fcr + 0:2^^) + ai ai r{k + r) 

+ {1 - e) (^(^- aiai k{k + £ + r) -al k{k + r) + ai a2 {k + £){k + £ + r) 
+ 02 Oil {k + £){k + r)) (fc + rf 

2{k + rf 



- ai ai (ai k{k + r) - as (fc + ^)(fc + r)) 
4 



+ 



- ai^ (a2 a£ (fc + £ + r)^ - a\ t) 



a\ a^Dx -D2 -D3 

+ OLx {k + r){k + £ + r) (^ae {k + £ + rf + ai f - ai ai j + (fc + r) 
- «! ^2 (fc + + ^^)^ + «! (fc + ^')^ + ai ai {k + r){k + £ + r) - 



(65) 



(66) 
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ai a\ ( ai Q"^) ] 

Oil ' 

+ q;i(1 -e) ax oi{lt - axai{k^r){k^ t^r) {oLi{k^ t^rf 

+ (/c + rf [{-2 ai ai + ^2) {k^i + rf - af {k + rf 
+ ai (-«! + ai- aa) {k + r){k + £ + r) + ai al ( Q^) 



"^"^14 



al{k + rf{-{k + Tf + ^-a,Q^) 

Oil 



(67) 



+ ^(l-e) (^aiai ([k + rf + a^Q'' ) ( 2(fc + + r) + r(fc + r) + ai Q^) 

+ (fc + r)' {ai^ (fe + + £ + r) + a? (fc + r)' + ai al Q") ^ 

+ ai ai (fc + r)2 - + £) + - r) (fc + £ + r) - r^^ 
+ ai (-(fc + rf - ai Q') (a^ ( 2(fc + £)(fc + r) + r(fc + r)) 

ai {2e-r){k + r)^ 



+ 



+ 

2(fe + r)2 



" , (ai^ + a^') (xi {k + r)2 

q;iQ;^L)iL)2 -D5 
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2(ai^ + Q;2^ )(fc + r) 



23 



+ k{k + i) {k + r){k + i + r] 



AA22 



-ai + ai)[ - k{k + r) {k + £){k + £ + r) 



(68) 
(69) 

(70) 



+ (a, {k + £){k + r) + ai{k + rf^ (^^-^^^ fe(fe + £ + r) 



a2 {k + £ + rf + 



a2^ai {k + rf £ 



+ 



ae / 



+ e 



CKl Oil CC^ 

ai k{k + £ + r) + a2{k + £){k + £ + r) 



(ai - a2) 
ai 



■{k + rf 



+ 2{k + r){k + £ + r) ) 
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+ - aik{k + r) + a2 {k + i){k + r)^ - cca (fe + ^ + rf 



2{k + r){k + £ + r) + aiai C + 



^ A(k + r){k + £ + r) ^^^^ 

0.2 0L\ 
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24 — -r^ 2 



(71) 



+ ai - ai + r) + ^2 (A; + + r)) (^(fc + r)^ + ax )^ 

+ ai (fe + r)' (ae k{k + £)-a2{k + if + ai k{£ - 2 r) + ^l^li! 

+ a2{2£-r){k + £)^ 
- a2 (^{k + rf + ai (a^ {k + £){k + r) + ai {k + rf + ai{£ + r) {k + r)) 



+ 
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aiaeD2 D4D5 



{k + £ + rf e 



33 



•34 - - A^c 



2 e ( — ai q;^ ((fc + r)(fc + ^ + r))' 

Q!2 



Z^l 1^2 /^3 

+ (fc + ^)(fc + r) (^-Q;2Qi^(fc + -^ + r')^ + Q^iQi^('^ + ^)^ + «i«i-^^) 

+ CKi a\ + (fc + r)(fc + ^ + r) fai (2 ai ^2 + 012 (-«! + ^2)) {k + r)^ 

ai 0:9 V 



(72) 
(73) 
(74) 



+ 2 cti q;2 £{k + ^) - «! q;2 a? 



+ «! (fc + r) 



(ai - eta) (fe + ^)(fc + £ + r) 



2 I --, „, /n2 



+ Q;£(fe + ^ + r)^ 



+ 



3 (ai - ai) ai (fc + rf r{k + £ + r) 



+ ((ai - ai) (fe + ^)(fc + r) + ""^^ + ai a, ^2 Q") (£ + r) (fc + ^ + r) 



^ — (^af ae{k + £ + rf — 0:2 cci (fe + rf — ai 0:2 cti £^ j 
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• (2 ai at {k + i){k + r) + ax at {k + rf + ai ax {£. - 2r) {k + rj^ 
- {k + £){k + £ + r) (ae {k + r f + ai 612 + (ai - ax) [i + r) (fe + r)) 

+ (fe + r)(fe + ^ + r) ( {a-i OLx - ax a^ a2 + 2a2 {ax - 0(2) ai 

\ax a2 

— 4 tti 0:2 ai) (k + rY — 2 ai 0:2 Oi£ 



+ 



{—ax + ax) 

{—ax + OLx ) 
a2 



{-at {k + if + axl^ - ax at Q^) 

-(2 at {k + £){k + r)- ^2 {i + r) {k + £) 



-2ax{e + r) (fc + r))) 



+ 



4{k + r){k + £ + r) ^^^^L 

a2 ax 
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Nr. 



■ ax {a2 ax + ax a2) {k + r){k + £ + r) 



aiDx D2 D3 

-(1 - e) (ax^ {k + £f {k + rf + a? (fc + rf 



2 


\2a^ 


D2^D^^ 


ax ax 



+ ax^ax {k + £){k + r) 

+ ax ax (2 {k + rf + axQ^) {{k + £){k + r) + ax Q^) 
+ ax[-ai{k + £f {k + rf + 2 ax {k + rf ((^ - rf + r^) 



+ ae{k + r)^ - 



ax {k + ry 
ax 



-axQ'' + 3{£-r) + ) 



+ 



3 ai (^{k + rf + ax Q') (ae {k + £){k + r) + ax {£ -r){k + r)) 
2(fe + r)2 



A4^ 



ax' 



44 



Nr 



(ai^ + a2^) cti (A; + r) 



ax D2' 



(ax k{k + r) + ax{k + r){k + £ + r)^ 



N, DxD2DeDr 
■ (^{k + £f-2ax {k + £){k + r) + a2axQ^') 

+ (^ax k{k + £) + ax{k + £){k + £ + r)^ (^-2ax{k + £){k + r) 
+ {k + rf + ax a2 Q^^ 
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+ e(^{k + £){k + r) (^{ai ai + as "2) + ^ + r) - 2 ai ^2 + r) 

- 2 tti 0:2 (fe + + ^ + r) + 2 q;i ai (fe + £){k + r)) 
+ (-«! ai + aa "2) + r) {k + i){k + £ + r) 

- k{k + £){k + r){k + £ + r)) 

+ ai{k + rf (k{k + £) + ai on Q^) 

+ ai{k + £f (^{k + r){k + £ + r) + aiai Q^^ 

+ ai ai (ai k{k + £) + ai {k + r){k + £ + r) + ai ai Q^) ^ 



L»i L»2 D-i L»7 



— (q;2 «! — 2 «! tti a;2 — «! "2) (fe + + £ + r) (fc + r) 



Oil 

(80) 



+ — (fc + £)(fe + r) fai (q;2 ai + ai ^2) ik + £ + rf 

ai V 

+ «! («! — 0:2 ) 0:2 (fc + r)^ + «! «! 0:2 (4 {k + r){k + £ + r)-a2Q^ )) 

^^/ Ka^ + 2c.,a,a2-C.2^i) (fc+^)(fc + ^ + ^)(fc + ^)2 

V Oil 
+ ai al{k + £){k + £ + r) 

+ —{k + £){k + r) fai (ai cti + q;2 ^2) (fe + ^ + rf 

- ai (cKi - 0:2 ) "2 (fc + r)^ + cci ai 0:2 (-4 (fc + r)(fc + £ + r) - cti Q^))^ 

(^ai fc^ - q;2 fe(fe + £ + rfj (a2 {k + £){k + r) - ai {k + rf^ 

(81) 

+ {ai^ - a2^) (k{k + r) {k + £){k + £ + r) - k{k + £) {k + r){k + £ + r)) 
-ai ( -k{k + £ + r) + {k + £f + 

■ [cxi {k + £)(A; + r) + {k + rf + ai ai Q^) 

+ e 1^ - ai - ai fe^ + q;2 fe(fe + £)) (k + rf 
+ (ai'-Qi2' )(^-k{k + r){k + £){k + £ + r) 



1 2 
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+ k{k + £ + r){k + £){k + r) - k{k + £){k + r){k + £ + r)) 

2 — 

+ (k + r)(k + £ + r) (-a2aie + ai^ (k + £f + ^^^-^ £''] 
ai e + at £{k + ((fe + r)^ + Q^^ 



+ 



2 -^^13 



H — -AB^ 



-2e(—a^ae[{k + £){k + r)]^ (82) 

+ ^^i^i:i^^^(/c + + € + r) (fc + r)^ 

-ai{k + r){k + £ + r) (^{k + £f + on Q^^ 

- tti ae (oLx {k + £f + ai (fc + r)^ + ai ai j 
+ q;^ (fc + + r) (^2 - r) (fc + ^ + r) 

+ + ^lt^i±^) (fe + - a? 1 

- ((fc + r)' + ai ^2 <5') ( - (fc + £)(fe + ^ + r) + ai - 2r) (fe + 

- — (fe + r)^ f - tti (fe + If - ^2 a? OLi + «! (-«! + tti ) + r) (fe + ^) ) 

+ (fe + + r) [ - ^i(-^i + ^i)(^2 + 2c.,) ^ 
L q;2 

+ ^^i^(^(fe + rf + a2Q^) 

+ f «2 (-a^ (fe + + ai ^2 _ g2 + _ 2r) {k + £ + r)) 

-2ae {ae {k + r){k + £ + r) - ai{£- 2r) {k + r))) 
+ (-ai + ai) (^{k + r){k + £ + r){£- 2r) (fe + £) 
- (k + £)(k + £ + r)(£- 2r) (k + r)) 

+ {{k + £f + a2 ai Q') (2 ae{k + r){k + £ + r) + ai {£ + r) {k + r)) 
a2 ai 

2 

- N, — ai (a2 ai + ^2) (fe + £)ik + r) (83) 

q;^ Di D2 D5 D-j 
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1 2 


(ai - ai ) 


Nc D2DsD4D7 





-a2 Oi\ + «! q;2 ) 

Qfl 



(fc + rf (a2 k{k + £) + a2{k + £){k + £ + r)) 

(84) 

{k + £f (a2 k{k + r) + a2{k + r){k + £ + r)) 
-2{k + £){k + r) (al k{k + r) + (fc + + + r) 
+ a2 6i2 {k{k + £) + {k + r){k + £ + r))) 

+ e(^{k + £){k + r) (^{ai ai + a2 ^2) k{k + £ + r) - 2ai a2 k{k + r) 

- 2 ai a2 (fe + + ^ + r)) +2 as ^2 [(fc + + r)]' 
+ (ai ai - a2 ^2) + r) (fc + + £ + r) - + £) {k + r){k + £ + r)^ 
+ (ai ai - a2 ^2) f — fc(fe + £)(k + ^ (fc + £)^ (k + rf 



+ 



— ik + £f ik + r)ik + £ + r)) 
a^ / 



ai D2 ^4' ^6 



Q;i^(l-e) aife^ -Q;2(fc + £)^ 



q;2 «! -^^ 



^(fc + r) (85) 



+ 



Oil 



Nr. 



1 



D2D4D5 Dr 



2e 



0(2 {—Oil Oil + CK2 Qi2) Q;^ 



(fc + ^)^(fc + r)^ 



(86) 



+ aek(k + £) (2 (k + £)(k + r) + ^^^^^^^ (k + r?) 

- ae k{k + r) (^{k + £f + 2{k + £){k + r) + ai on Q^^ 

- (k + £)(k + r) ((k + £f-2(k + £)(k + r)-ai on qA 
+ ik + £)ik + r) (^»2»dik + £r + a2a,Q') 



+ 



(«! -ai) 



{k + £f {k + rf e 
0L2 \ z \ Q ) 

0L\ OL\ OLi 



+ ^(-2 at k(k + £) + aik(£-2r) + 2 a2 U - 2r) (k + £))) ] 
+ {oii - ai) - k{k + r){£- 2r) {k + £) + k{k + £){£- 2r) {k + r) 
-—{k + rf (2aek{k + £) + a2 {£ + r) {k + £)) 
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+ J_(fe + £)2 (aek{k + r)-a2 - 2r)(fc + r))^ 
+ ^ ((fc + if + a2 «i g') {k + rf + ai + r) (fe + r)) 

Q;2 oti 
2 

Ar<3 — ^ ^ ^ ^ a2 {a2 ai + ai ^2) (fe + + r) 

at D2 V4 D-j 



2 



AT, DiDsD^De 



(ai {k + £){k + £ + r) + ai{k + r){k + £ + r)) 



■ - fc^ + 2 cci + £ + r) + CKi 0:2 ( 
(ai k{k + £) + ai k{k + r) j (2 k{k + £ + r) - ai{k + £ + r) 

,(fc + r)2 e 



+ q;2 ai [- — ■ — — H )] 



+ e (^(ai ai + a2 012) k{k + £ + r) {k + £){k + r) 

+ 0:2 - 2 k{k + £ + r) k{k + r) + k"^ {k + r){k + £ + r)^ 

+ a2 a2 (k{k + r){k + £){k + £ + r) - k{k + £) {k + r){k + £ + rj^ 

+ aia2(^-2k{k + £ + r){k + £){k + £ + r) + k{k + £){k + £ + rf^ 

- ai ai (k{k + r) {k + £){k + £ + r) - k{k + £) {k + r){k + £ + rj^ 

+ ^f-a,ik + r)ik + £ + r) + a,e) + ^) 

- ai ai - 2 {k{k + £ + r)f + k'^{k + £ + rf 



Nr. 



2 e (^a^ (^ai (fe + - ^2 (A; + £)(fe + r)^{k + £ + rf 



{k + £){k + £ + r) {k + rf 



-ai{k + r){k + £ + r) ({k + £f + alQ'' ) 

+ ai ae (^^^^ + ^) («i {k + £){k + r) + a,{{k + £f + a, Q^)'^) 
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_ a, ('I^±ll! + ^) (a, a, {k + + r) 
V Oil a(,/ \ 

+ 2axOLt [{k + if + ax Q^) + ax ai (£ - 2r) (fe + 
+ (fe + ^ + r)2 - 2 ai (fe + if + as (fc + + r) 

+ ai a2((^-2r) + 
+ 3 (ai - aa) ( - (fc + r)(fc + £ + r) r(fe + £) + (fe + £)(fc + r) r(fc + £ + r)) 
+ ((fc + if + Q') (a^ (fe + + £ + r) + ai (£ + r) (fe + £ + r)) 

+ (fc + £)(fc + £ + r) { 2axai^{k + rf + ax^axi^ 

a\ a\ ai \ 

- ax ai [{k + if + OLx - a^ a^ a^ ({k + £)(fc + r) + i{k + ^) )^ 

I 4(fc + £)(fc + £ + r) ^ 

ABl = TV, ^ «i (a? + al) (k + i){k + i + r) r' (90) 

ai Di D3 D-j 

- -4af{{k + i){k + r)f+^-^^^^^^^^^{k + ifik + rf 

(91) 

+ 4 e ("4 a/ ((fe + £) (fc + r)f - ("^^ ^ ^2) V ^ ^ ^^2 (93) 

+ a/ (ai {k + £)^ + ai{k + rf + ai ai Q^) 

(fe + £)(fc + r) (^ai {-ai + ^2) (A; + if + ai (-cti + ^2 ) (fe + rf 



abL = - a^c — 



ai ai 



- ai ai (cKi + CKi) ^ 
+ ^^{k + rf(a2aiQ^ - 3 {i - r) {k + i)) 
+ ^(fc + £)2fa2Q;,g' + 3 (i - r) (k + r) ) 



+ (fe + £)(fc + r) (^-2^ac^2ai + a,a2 + 2(a,-a2)a,) 

V «! 

_ 2 a^ (a2 «! + «! 0:2 + 2 (ai - 0:2) a^) _^ ^^2 

+ 4 (a2 + ^2) af + 2 (aa ai + ai ^2) ( - + 2 + 2 (£ - rf^ 
- 6 (ai - a2) ae {i-r){k + i) + 6 {on - 0:2) ae {i -r){k + r) 
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+ 3 ai Q'^ {£ - r) (k + r) - 3 ai 612 {I - r) {k + i) 

ABl^ = N, («2 «i + «i «2) (fc + + r) (93) 

^^^5 = (94) 

At several places we have expressed the matrix elements AAj-j-^, ABj'^js^ in terms of the 
corresponding matrix elements for the longitudinally polarized photon; we therefore list here 
also the expressions for longitudinal polarizations 0: 

2 

AA^^ = Cf -f:;^ (1 - e) ai ai Aj (95) 



DjD 



2 



AA^2 = ^ D njD "2 ^ "2 "1 '^^) (96) 

AA^, = Cf j^^l^j^^ ((l-e)aMi^3 +«? (fc + r) {k + l + r)) (97) 

^^^4 = - A^c ^2 "1 ("2 ai + 2 (1 - e) afj Ai A4 

1 1 \ 

H — aia2Q;i-D2H ai 0:2 ai Air — ai 0^2 Q^i ^ (98) 

2 2ai J 

-^22 - Cf _ n; n2 (1 - ^) "i °2 ^2 (lOO) 

«! 1^2 

•^2^ = - ITT n n n ((^ ~ ^) '^2 "i ^2 ^3 + ^2 "i (A; + r - «! r) (fc + £ + ai r)) 

I\c L>1 L>2 ^3 -t^4 

(101) 

AA^ = Nc ^2 n ( o ^ - "1 "2 («2 ai + 2 (1 - e) a^) A2 A4 

11 \ 

H — aia2CiiD2A ai ^2 'Si ^2^" + «i «2 'Si ^ (102) 

2 2a^ / 

= - AT, ^2^n n "1 "2 tti (103) 
ai D4 



AA34 = Nc ^ ^ (77—^1 ("1 «i - «£ - 2 (1 - e) a|) A3 A4 

— - ai ^2 -D2 + oi\ a2 + - — ai 0^2 tti Aar + ai 02 r 

2 2a£ 

+ ^a2ai A3(fc + £) j (105) 



(104) 



22 



AAl, = - TVc n n n n «2 (106) 

L>i L>2 L>3 JJ5 

AA^ — Nc ( — A4 «! ((1 — e)a'j + a2 cti) + ai aiD2 + 2 cti al cti r^] 

(107) 

AA^^ = - A^c «2 ai (108) 

^^55 = (109) 
A-Bfi = ^ ^"^^ ^ ((1 - e) ai ai Ai Si + a, cxi (k + aiv) {k + i + as r)) (110) 

AB^2- -Cf ^ ^ ((l-e)aia2AiB2 + aia2 (fe + r) (fc + ^ + r)) (111) 

iJl iJ2 -/^S i->'7 

= K d^dUd, ( - ^' - ("2' 

1 /I 

Ai3f4 = -— ai (tti ai - - (1 - e) 2 aj) Ai 

+ ai ai a2 D2 — —cki 0:2 L/7 + - — ai 0:2 cci q;2 AiT* — ai 0:2 cci q;2 t* 
2 2a£ 

-^aia2Ai{k + r)^ (113) 

ABfg = - TV, ]^ a, a2 a\ ^2 (114) 

ai Di D2 Dr, D-j 

■^^22 = TT ((1 ~ ^) 0:2 A2 B2 + 0:2 0:2 (fe + r - ai r) (fc + £ + cti r)) 

l\c D2 J^3 J-' A J^7 

(115) 

AB^, ^ - Cf _ ^ (1 - 6) ai a2 A2B, (116) 
ai V2 Ui JJe 

AB^^ = - Nc ^ ( TT- «2 (ai ai - - 2 (1 - e) a^^) A2 B4 

+ ai a2 Oil C(2 D2 Oil d-l «! -D7 H cti 0:2 cii 0:2 Asr + ai 0:2 cSi 0x2 f"^ 

2 2 q;^ 

- ^Q;2a2A2(fe + r) ^ (117) 

AB^^ = ^ ^ ai al Oil 0i2 (118) 

D2 D4 L>5 L>7 

^1^33 = TT n n n n ((^ ~ ^) ^3 ^3 + (fc + r - Cti r) (fe + £ + r)) 

(119) 

■^"^ = - aa'a,c, (li^"' («■ a. + 2 (1 - .) .a A3 

H — cci a2 D'j H cei a\ 0.2 A2,r + cci 6l\ 0.2 r\ (120) 

2 2ae, J 
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ABg = a, al al (121) 

ai IJi L>3 U-j 

ABl^ = -N, ^ ^ ( (2 ai ai + - 2 (1 - e) a^^) A4 
L>2 L>^ Uj V 

1 1 

+ -ai ai (X\ 6l2 D2 + -ai ai ai Dj + 2 ai ol\ — a'j (k + £) {k + r) 

(122) 

AB^^ = Nc n2 n "2 tti "2 (123) 

JJ2 L)^ U7 

AB^, = 0. (124) 

Here we have used the following abbreviations: 

Ai=a2i + ai{k + r) (125) 

A2=aie + ae{k + i) (126) 

A3=aie- ae{k + r) (127) 

A^=a2{£-r) + ae{k + r) (128) 

Bi=a2£- ae{k + i) (129) 

^2 =ai£-a£(fc + r) (130) 

S3 =ai£ + a^(fc + £) (131) 

Bi=a2ii-r) - ai{k + i) . (132) 

We conclude our listing of the squared matrix elements with a few remarks. First, similar 
to the longitudinal case, the sum of the squared matrix elements for the transverse photon 
is ultraviolet finite: individual terms diverge as £ becomes large (for fixed k), but in the sum 
these divergences cancel. Secondly, it is easy to show that the factorization properties in 
configuration space which for the longitudinal photon have been discussed in section V of 
also hold for the transverse photon: this supports the validity of the photon wave function 
picture beyond the LO approximation. Finally, the limit of small ai agrees with the LO 
BFKL calculation; we shall return to this question further below. 



V. SINGULAR LIMITS OF THE REAL CORRECTIONS 

In order to formulate finite expressions for the qqg contribution to the photon impact factor 
we examine the divergent limits of the squared matrix element |A^gggP for the process 
7* + g — ^> qqg + q. For the phase space of the gg^f-final state we have from (0), (|l^): 

sdjSr da d^^'^k dae d^~'^£ 

(i<Pm^ = 2«(l -a-ae) (27r)^-i 2^ (2^0^ ^ ^ 

In the following we suppress the superscript {D — 2) in the integration measure of the 
transverse momenta. The contribution of the qqg state to the impact factor reads: 

da f dk f dae f di (o) ,2 . ,^ 



I 



2a{l -a-ae) J (27r)^-i J 2ae J (27r)^-i 
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where the square of the vertex function, |r^?_^qggp, follows from ( |62D and (|63|) . 

We begin by examining the limit in which the gluon is either collinear to the quark or 
to the antiquark of the qqg-sjstem. In the latter case, the gluon momentum i is collinear 
with the momentum q — k — i of the antiquark, resulting in D4 = or 

£' = £+-^k = 0. (135) 
1 — a 

We approximate our results for iFy^'^I^^gg^p in the region around the singularity, i.e. we cal- 
culate the residue of the pole at D4 = 0. The collinear limit of QCD matrix elements is well 
known and can be taken from the literature . As an additional check, we have calculated 
the divergent limits directly from our matrix elements. For both polarizations of the photon, 
the g \\ q collinear limit of (|134[ ) can be written as 



27r 



p{o) ^2 da dk 



in co\\,q 
.(0) 1 2^^ "5"Pr 

\ dai d£ 



2a(l - a) (27r)3-2^ 



X CpPna — e — ^ ^ (136) 

Ua,k)dadkC,P,, (y^,^) (2^)fi^.2 (137) 



where X2(a, k) = X2(a, k; r, Q) was defined in (^), (p8|), (p9|) , and CpPgqiz, e) denotes the 



q ^ g Altarelli-Parisi splitting function in D = 4 — 2e dimensions: 

. X 1 + (1-^)^ , , 

Pg,{z, e) = ^ - ez . (138) 

The color coefficient Cp indicates that, for the examination of singular limits it will be 
convenient to use the color basis Cp and Ca rather than the color factors Nc, Cp, and 1/A^c 
which appear in the matrix elements of the previous section. In this basis, the Ca part is 
free from collinear singularities. 

In the case where the gluon is collinear to the quark, the infrared limit is obtained from 
( |137|) by substituting^ abya; = l — a — and k hj k = —k — r — £. The "soft-collinear" 
limit of (|137| ) that is needed later is obtained by retaining from the splitting function only 
the pole 2/z: 



I (0) ,2^ , _ ^ 

\i- y*^qgg\ ^^Jqqg 



soft,coll,g , ,^ 

= Ua,k)dadk2Cp^ (139) 



Next, we examine the 'soft limit' where the gluon momentum — > and £ ^ 0. 
Formally we substitute £ ^ pi or £ ^ p£, pai, expand around p = and keep only the 
leading term. As in the collinear case, the result is proportional to the Born approximation. 



^This is nothing but the (? <-> g variables exchange at fixed gluon momentum. 
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for both polarizations of the photon. We separate the result according to different colour 
structures: 



|p(0) 1 2 



soft 



Ca ir 



(0) ,2 



soft 



Ca 



+ Cf ir 



(0) ,2 



soft 
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Ca ' (£-a,a)'(£ + a,b)' 



_ 1^(0) 12 4a^^(a + b)- 



where, for convenience, we have defined 

k + r 



{£ - apay {£ + aeb) 



k 



2 ' 
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(141) 



(142) 



a 



a 



1 — a 



(143) 



In the soft limit, the phase space measure (|13 
(p5|) and differentials of the gluon variables 



factorizes in the product of the qq measure 



d(j), 



sdPr 



119 



271 



sdPr 



dai di 
2ai (27r)^-i 



(144) 



Our result for the soft limit agrees with [^]. We also note that, in the limit i' = £ + aeb 0, 
we reproduce from ( |142| ) and (|144|) the "soft-collinear" limit ( |139| ). 

Finally, we have to consider the limit in which the gluon is emitted in the central rapidity 
region. Formally, this corresponds to the singular limit 0: the matrix elements listed 

in the previous section are finite for ae —>■ 0, but the phase space element (^^ introduces 
a logarithmic singularity. Consequently, when doing the subsequent ae integration this 
divergence ~ 1/a^ will produce the leading-Ins logarithm. As discussed in section |T|, this 
part of the gluon phase space belongs to final states with two large rapidity gaps and has 
to be subtracted from our result. This will be the topic of section For the moment we 
only state that in the limit — 0, for both polarizations of the photon, the sums of our 
matrix elements split into the product of the Born approximation and the BFKL kernel: 



-(0) I 



<0) 



LL 



*qq\ 



I ACat^ {(- - r) 



(145) 



Here the last factor is the BFKL kernel (without the gluon trajectory function) 

Ca (£ - rf 



}CreJr-i,r) 



TT 



£2 



(146) 



and 1/ (£ — r) are the propagators of the (reggeized) gluons. Note that in the Born level 
vertex functions the transverse momenta are shifted from r to r — £ and from k to k + £. 
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The color factor Ca signals that the Cp part of our matrix elements does not contribute to 
our subtraction of the central region. In the limit where the centrally emitted gluon is soft 
{i <ti r), we obtain 



|p(0) I 
I J- ^*^qqg\ 



soft 
LL 



.(0) 



'a,k)\ 



£2 



(147) 



The same result is obtained from the soft limit (|141|) , if we go into the central region limit 
-C 1. 



VI. SUBTRACTION OF THE CENTRAL REGION 

After having calculated the various divergent limits of our matrix elements, we now turn to 
the subtraction of the central region. As we have discussed in section ||, contributions to the 
total cross section are divided according to the number of large rapidity gaps. What we have 
presented in the last two sections are the results for the process 7* + g — (qqg) +q, imposing 
the constraint of a large rapidity gap between the [qqg) system and the quark. These results 
therefore contain, as a special case, still the configuration where the qqg system contains a 
second rapidity gap between the qq pair and the gluon: this piece has to be removed since it 
counts as a configuration with two gaps and, in order of the total cross section, it entirely 
belongs to the leading-ln s approximation (|1^) . 

In order to make the separation more explicit, we first divide the phase space of the 
produced gluon into 'upper' and 'lower' halves. Introducing 

at' = ^, (148) 

the 'upper' region a™* < ai < 1 — a consists of the 'upper' half of the central region plus the 
fragmentation region of the incoming photon, whereas the 'lower' region ae < a^"* contains 
the lower half of the central region plus the fragmentation region of the quark. In the upper 
region we write: 



I ^\*^qqg I d(j)qqg 9 (o^ - O 

>1 



dkd£ da dail3{ae,£;a,k;r,Q)Q{ae - a'^"') (149) 
Jo Jo 

Here X3 denotes the full integrand involving the sum of the matrix elements of the three- 
particle final state. 

The division into central region and fragmentation region is closely related to the choice 
of the energy scale sq. In (ITSD we have introduced the energy scale sq: since a change in 



So can be absorbed into the NLO impact factors and one might expect that there 
exists some freedom in choosing this scale. We adopt the simplest choice, namely a constant 
value which is independent of kinematic variables inside the impact factors. Examples for 
So include the (negative) virtuality of the t-channel gluon, r^, or the (negative) mass of the 
photon, Q"^. 
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In order to define the qqg contribution to the photon impact factor, , we have to 

subtract, starting from ( |149D , the leading-Ins piece with the energy scale Sq. This can be 
achieved by considering the central region limit (|145| ) of the matrix element which implies 



X3(a£,£;a, fc;r,Q)|LL =X2(a,fc + £;r - l,Q)- 



2tx 



\2-2t 



(£-r) 



(150) 



We recognize the integrand Xi occurring in the LO impact factor (^) with the value of 
the reggeon momentum being shifted to r ^ r — i (since the momentum t is carried by 
the outgoing gluon), the propagators associated to the t-channel (reggeized) gluons, the 
measure factor l/(27r)^~^^ required by the definition (|]), the real part of the BFKL kernel, 
and, finally, a factor 1/a^. It is this factor 1/a^ which, when the integration over is done, 
provides the logarithm of the energy. Our aim of finding the LL piece with energy scale Sq 
leads to the following ai integral: 



e - aT) = - log 



(151) 



with af^^ = \i\/^/so. Therefore, we define as genuine real emission contribution to the 
photon impact factor the expression 



(l.real) 



dkd£ da da£ 

Jo Jaf-^ 

l2{a,k + i;r -l,Q) 



Isiae, i; a, k, r, (5)6(1 - a - ai) 



]Cro.Jr-Lr)—e( 



(27r)2-2e(^ _ ry 



■ (152) 



In this definition of our impact factor, the following remark is in place. With our choice of 
a constant energy scale, which has led to the subtraction, we are adding (in the LL term) and 
subtracting (in $^l''^*^'^''*) a region of the phase space which is kinematically forbidden. In fact, 
the LL subtraction (last line of Eq. ( 152 )) includes the region of phase space where a^ > 1—a. 
It is important to make sure that this "forbidden" region \£\ > (1 — a)^yso , 1 — a < ai < 
\£\/-,/so contributes a finite term and does not introduce new divergences. The finiteness 

results from the "good" UV behaviour of the LO impact factor, namely ^^}{r —£) ^ const 
for £ — s> oo (up to logarithms). After multiplication with the propagator ~ l/£^ and with 
the factor log[|€|/(l — a)y^] (which is due to the ae integral), the resulting £ integral is 
finite. This situation is quite different from the case of partonic impact factors 0,0, where 
$g°g(r— €) ~ and the £ integral, after including the gluon propagators and the logarithmic 
factor, leads to a double pole in 1/e. 



VII. FINITE COMBINATIONS 

In the previous section we have suggested a suitable form of the LL subtraction and thus 
specified the real emission contribution to the impact factor. What remains are the infrared 
divergences which have to cancel once we combine real and virtual corrections. In this final 
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section of our paper we show explicitly the cancellation of all the IR singularities, and we 
present expressions consisting of well defined finite terms. 

Given that all the singular terms of the virtual corrections are proportional to the LO 
impact factor (^71) , as one can see from Eq. (0), it is convenient to study the singular 
regions of the phase space integrals of the real corrections by integrating first only over the 
gluon variables. In fact, the gluon integration generates singular terms which are in turn 
proportional to the LO impact factor. They will be combined with the singular virtual 
corrections after a proper identification of the remaining integration variables. 

It is also convenient to separate the different colour structure contributions into a Ca 
term and a, Cp term 

J3 = CaI^^ + Cpl^^ , (153) 

because the LL subtraction is contained only in the former term, whilst the coUinear singu- 
larities are present only in the latter. Soft singularities are found in both terms. 



A. Ca term 



Let us consider the Ca contribution to the real emission impact factor corrections, as spec- 
ified in Eq. ( |152|) . We want to compute first the I and integrals at fixed k and a: 



da dk 



/POO 
d£ / dai[I^'^{ae,£;a,k]r,Q)Q{l — a — ae) 



-l2{a,k + £;r -l,Q) 



(27r)3-2e(r-£)2£2^^ 



(154) 



This integral is divergent in the soft region — > 0, £ ^ when e = 0. We evaluate 
it in dimensional regularization by means of the subtraction method: we subtract from the 
integrand its soft approximation, in such a way that the resulting integral will be e-finite. 
Subsequently we re-add what we have subtracted, and we integrate this term analytically. 
Of course, the soft limit is not uniquely defined outside the soft region. We choose the soft 
approximation according to Eqs. ( p.41| , |147| , |144| ) ; we extend the soft subtraction up to the 
kinematic limit of the gluon longitudinal momentum, i.e., ae < 1 — a: 



d<^. 
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soft 
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(155) 
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Accordingly, we define the finite Ca part of the real corrections to the impact factor as 



^(l,real) 

7* 



finite 
Ca 



dkdi 



X2ia,k + £;r -l,Q) 



da 
2 



daA [I^^{ae,£;a,k;r,Q)e{l 
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1^1 



(156) 



6(1 



a 



at) 



Note that we have shifted to zero the lower limit of integration of ai because, thanks to the 
LL subtraction, there is no contribution to the integral coming from the central region. Note 
also that the singularities of the integrand due to the denominators i — a^a and I + aib are 
integrable, because of the scalar product in the numerator. This reflects the fact that the 
Ca part is free of collinear singularities. 

The expression written in Eq. ( |155|) can be computed analytically. Since the lower limit of 
integration in depends on the energy s, one gets an s dependent result. However, because 
of the LL-soft subtraction to the soft term, the central region contribution is suppressed, 
and we expect a finite expression in the s ^ oo limit. We can compute the integrals of the 
two terms in Eq. (|155|) in this limit (a^"* — > 0) separately because dimensional regularization 
promotes the divergent log s contributions into UV e-poles, which will cancel in the difference. 
We have: 
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0( 



TT 
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i-2e 
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These two integrals have to be multiplied by X2 and integrated over a and k. After Eq. ( P?]) 
we have noticed that I2 is symmetric under the change of variables a ^ 1 — a, k ^ —£ — k, 
corresponding to the q q exchange; the term (a+6)^ = M^/ a{l—a) is also invariant under 
this exchange. Therefore, on the right-hand-side of Eqs. ( |157| ) and ( p.58| ), we can replace 
(1 — a)~^^ \(^~'^^ + |(1 ~ ct)"^*^) which casts the expressions into explicitly symmetric 
forms w.r.t. a ^ 1 — a. This manipulation will be useful when combining real and virtual 
corrections. 

In conclusion, the divergent Ca part of the real corrections to the impact factor is 



7* 



(l.real) 



divergent 
Ca 



= / dk 



daX2{a, k] r, Q) 



(27r)3-2^ 
cr 



{Ji - J2 



J dk J da l2{a,k]r,Q)j^^i^^ [log M'^ -log So -log a{l - a)] 



(47r)2-^ 

+ [(logSo + loga(l-a))'-log2M2] +0(e)| 



(159) 



For consistency, the variables a and k of the above equation have to be identified with the 
analogous variables of Eq. (^) . 
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B. Cp term 



Let us now consider the Cp contribution to the real emission impact factor corrections. The 
first thing to note is that, since the BFKL kernel is proportional to Ca, the LL subtraction 
of Eq. (|152|) does not affect the Cp part. Therefore, the whole Cp part of T3 contributes to 
the impact factor. In addition, we can shift the lower limit a™* for the ae integral to zero, 
just because the central region does not contribute in this case. 

The Cp term contains both soft and collinear divergences. We separate the singular con- 
tributions by means of the subtraction method. In what follows, the last two parameters of X2 
are always r, Q. Therefore we simplify the notation by indicating T2{a, k) = X2(a, k; r, Q). 

We start by identifying the soft singular term. This is done by subtracting from if^ its 
soft limit and then re-adding it again. We remark again that, beyond the soft region, the 



limit of integration is not uniquely defined. In practice, by using the expression ( |142| ) , we 
adopt an asymmetric function with respect to the q q exchange at fixed non-vanishing 
gluon momentum: 
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where in the last equality a symmetrization in a, as indicated after Eq. ( |158| ), has been 
performed. 

Both the full Cp term oc X^^ and its soft limit (^6^) have a divergent behaviour when 
the gluon is either collinear to the outgoing quark or antiquark. Also for the evaluation of 
the collinear divergences we adopt the subtraction method. For the case of the gluon being 
collinear to the antiquark [g \\ q), we derive the limiting expression of the original integrand 
(i.e., without soft subtraction yet) from Eq. ( |137D , which yields {£' = £ + agb) 
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(161) 



where we have decided to perform the collinear subtraction only in a cone — specified by 
the G function — whose axis coincides with the gluon momentum being parallel to the 
antiquark momentum, and the vertex sits on the soft point. The "opening angle" of the 
cone is parametrized by a cutoff A: the bigger A, the wider the domain of the collinear 
subtraction. 

In order to end up with a finite integral, we have to subtract also the collinear singularities 
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of the soft subtraction term ( |160|) . For the g \\ q we have (see Eq. ( |139| ): 
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The soft-colhnear subtraction uses the same cone as defined for the colhnear one. It is 
important that Eq. ( |162| ) can be obtained both as a colhnear hmit of Eq. ( |160|) and as a 



soft hmit of Eq. ( |161| ). This guarantees not only that the colhnear singularities introduced 
by the soft subtraction are subtracted, but also that the soft singularities introduced by the 
colhnear subtraction are properly removed. 

Note that the colhnear term ( p.61| ) is proportional to X2{a,k). This is a consequence of 
the convenient choice of variable employed to parametrize the outgoing particles. In fact, 
when the gluon is colhnear to the antiquark, the gq system has to be treated as a single 
entity carrying the sum of quantum numbers and momenta of the constituting particles. In 
this case, the gq longitudinal momentum fraction and transverse momentum are 1 — a and 
—k respectively. The integration over the gluon variables can therefore be independently 
performed leaving the LO integrand I2 factored out. 

The treatment of the g \\ q colhnear singularity can be performed in the same way as for 
the 5' II ^ colhnear case. There are only two small differences: first of all, the colhnear limit 
of the integrand X3 is proportional to X2(a + a^, + €), which involves explicitly the gluon 
variables and I. This is due to our set of variables that does not describe the qg system 
and the single gluon independently. A more suitable set of variable in this case is simply 
obtained by changing the q ^ q labels: 
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(163c) 



From Eq. ( p.61| ) one gets 
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The additional factors in the second line on the rhs of Eq. ( |164| ) (when compared to the 
corresponding line of Eq. (|161|) ) are due to the fact that 
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(165) 
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The second difference is that, considering the soft-colhnear subtraction, we cannot use 
the soft hmit of Eq. ( p.64| ) at fixed a and fc, because this does not correspond to the colhnear 
hmit of the soft approximation ( |16CI| ). The right way of subtracting the soft singularity 
of ( |164| ) and at the same time the g \\ q colhnear singularity of (|16CI|) is to take the colhnear 
limit of the latter: 
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In conclusion, we write the finite part of the Cp term of the impact factor by subtracting, 
from X^^, the unintegrated expressions of Eqs. (|16CI| , |161| , |162| , p.64| , |166|) : 
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The divergent part of the Cp term is obtained by re-adding the subtracted pieces in their 
integrated form. In doing this, we identify the 1 — a and —r — k variables of Eq. ( |164| ) with 
the variables a and k of Eqs. ( |161| , [L62| , |166| ) (and later with a and k of Eq. (|3^)), so that 
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The final expression for the one-loop correction to the photon impact factor can be 
obtained by summing the real contributions given in Eqs. ( |156| , |1 59| , |167| , |168| ) to the virtual 
corrections $^1'^^'^*'^^''' of One can check the cancellation of soft and colhnear singularities 
by combining the singular piece of the virtual corrections (^) with the e-poles of Eqs. (|159|) 
and (|168|) . The cancellation of double and single pole is straightforward. 

In conclusion, the final expression of the impact factor at NLO can be obtained by 
summing the finite part of the virtual corrections, the scale dependent term (44) and the 
extra term ( |^ ) as discussed in section |T^, the finite terms of Eqs. ( |159| ) and ( |168| ), and the 
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integrals in Eqs. ( |156|) and ( |167| ) (evaluated at e = 0): 
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VIII. CONCLUSIONS 



In this third part of our program of calculating the NLO corrections to the photon impact 
factor we first have presented a list of the real corrections, the qqg intermediate state in- 
side the impact factor. This completes the results of the paper p| which contained the 
longitudinal photon only. We then have identified and computed the divergent parts of the 
phase space integral of the real corrections: in addition to the infrared singularities which 
are due to coUinear and soft configurations of the produced gluon there is a logarithmic 
divergence related to the Ins-piece of the central region. Removal of the latter part of the 
real corrections introduces the energy scale sq which represents a central element of the NLO 
corrections. Finally, we have combined the real corrections with the divergent parts of the 
virtual corrections: these infrared finite combinations are the main results of this paper. 

With these calculations we have essentially completed the analytic part of our program. 
What remains are numerical steps: the phase space integrals both for the qq state (which 
includes the finite pieces of the virtual corrections) and the qqg state. Certain parts of these 
integrals can be done analytically, but, as a general strategy, we define standard integrals 
and express our integrands in terms of these expressions. The final evaluations have to be 
done with the computer. 

Once this is solved technically, the photon impact factor allows us to compute the total 
cross section for 7*7*-scattering to next-to-leading logarithmic accuracy in s. This will be 
the first consistent calculation of a cross section to this order since the NLO corrections 
to the BFKL kernel were first published Taking advantage of the recently obtained 

impact factor for forward jets [|l^] in the same limit we may apply the photon impact factor 



also to the calculation of the forward jet cross section at HERA. 

We have emphasized already before, that interest in the NLO calculation of the photon 
impact factor comes from various directions. A prominent example is the photon wave 
function picture. In our analysis of the real corrections we have extracted the new qqg 
Fock component of the photon wave function which, most conveniently, is expressed in 
configuration space. In the present paper we have studied the cancellation of divergences, 
but this was done on momentum space. It will therefore be necessary to translate these 
calculations into configuration space. Furthermore, in order to obtain a clearer geometrical 
picture of how color charge is distributed inside the virtual photon, we need to take a closer 
look into the various pieces found in |0| and in this paper. Another question of interest 
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of the photon wave function picture is the form of the NLO correction to the qq Fock 
component and of the qq dipole cross section: do they destroy some of the main features 
of the LO calculations, for example the conservation of the transverse dipole size during 
the interaction with the target? Answers to these questions are to be found in our virtual 
corrections, but they require further theoretical efforts. 
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